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Abstract — Adaptive homodyne estimation of a continuously 
evolving optical phase using time-symmetric quantum smooth- 
ing has been demonstrated experimentally to provide superior 
accuracy in the phase estimate compared to adaptive or non- 
adaptive estimation using filtering alone. Here, we illustrate 
how the mean-square error in the adaptive phase estimate 
may be further reduced below the standard quantum limit 
for the stochastic noise process considered by using a Rauch- 
Tung-Striebel smoother as the estimator, alongwith an optimal 
Kalman filter in the feedback loop. Further, the estimation 
using smoothing can be made robust to uncertainties in the 
underlying parameters of the noise process modulating the 
system phase to be estimated. This has been done using a 
robust fixed-interval smoother designed for uncertain systems 
satisfying a certain integral quadratic constraint. 

I. INTRODUCTION 

Quantum parameter estimation (QPE) [1] involves estimat- 
ing an unknown classical parameter of a quantum system and 
plays important role in various fields such as quantum com- 
putation [2], quantum key distribution [3] and gravitational 
wave interferometty [4]. A common and technologically 
relevant example of QPE is estimating an optical phase. The 
fundamental limit to the precision of the phase estimate is 
set by Heisenberg's uncertainty principle [5]. On the other 
hand, the standard quantum limit (SQL) is the minimum 
level of quantum noise that can be obtained using standard 
approaches not involving real-time feedback. 

Since the phase of an electromagnetic field cannot be 
measured directly, all phase-measurement schemes employ 
measurement of some other quantity, that necessarily intro- 
duces uncertainty in the phase estimate. The standard method 
of measuring the phase of a signal is the heterodyne scheme, 
where the signal is combined with a strong local-oscillator 
(LO) field detuned from the signal, resulting in an introduced 
excess uncertainty scaling as 1/N (JV := \a\ 2 is the mean 
photon number). By contrast, homodyne scheme introduces 
greatly reduced uncertainty, in cases where there is some a- 
priori knowledge about the phase, by using an LO phase that 
is 7r/2 out of phase with the signal. Moreover, by adapting 
the LO phase using feedback during the measurement, it is 
possible to further reduce the excess uncertainty to obtain a 
mean-square estimation error lower than the SQL [6], [7], 
[8], [9] and even attain the theoretical limit [10]. 

However, these works were based on single-shot mea- 
surements of fixed unknown phase. Practically, it is more 
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relevant to be able to keep track of a time-varying phase 
instead [11], [12]. There are a number of ways of estimating 
a classical process dynamically coupled to a quantum system 
under continuous measurement, viz. prediction or filtering, 
smoothing and retrodiction [13]. Smoothing, in particular, is 
an estimation technique, that uses both past and future mea- 
surements, and, therefore, yields a more accurate estimate as 
compared to only filtering, that uses only past measurements. 
However, it is essentially a non-causal method that cannot 
be used in real-time but is used for offline data processing 
or with a delay with respect to the estimation time. 

The fixed-interval smoothing problem [14], [15], [16] 
involves measurements over a given fixed time-interval T. 
One solution to the fixed-interval smoothing problem is the 
Mayne-Fraser two-filter smoother [17], [18], [19], that uses, 
in addition to a forward-time Kalman filter, a backward-time 
Kalman filter, also known as an "information filter" [20] 
and finally combines the two estimates to yield the optimal 
smoothed estimate. Rauch, Tung and Striebel combined the 
information filter and the smoother into a single backward 
smoother [21], [22]. 

The first experimental demonstration of adaptive quantum 
phase estimation of a continuously varying phase using 
quantum smoothing was presented in Ref. [23], where an 
estimate could be obtained with a mean-square error of up 
to 2.24 ± 0.14 times smaller than the SQL. The experiment 
used a classical stochastic Ornstein-Uhlenbeck (OU) noise 
process to modulate the signal phase to be estimated. The 
authors have previously shown in Ref. [24] that the feedback 
filter used in Ref. [23] is only optimal when the noise process 
is a Wiener process and the measurement is assumed to 
be linear and that using an optimal Kalman filter, instead, 
significantly reduces the mean-square error in the phase 
estimate. Here, we show that using a Rauch-Tung-Striebel 
(RTS) smoother, in addition to the Kalman filter in the 
feedback loop, improves the accuracy of the phase estimate 
as compared to the (offline) estimator used in Ref. [23]. 

It is desirable to make the estimation process robust 
to uncertainties in the underlying parameters of the noise 
process, since it is physically unreasonable to specify these 
parameters accurately. Significant effort has been put on 
developing robust approaches to QPE [25], [26]. The authors 
have illustrated in Ref. [24] how the feedback filter and, 
therefore, the precision of the phase estimate can be made 
robust to uncertainty in one of the underlying parameters, 
based on a guaranteed cost robust filtering approach [27]. 
Here, the optimal RTS smoother is made robust to parameter 
uncertainty by applying the fixed-interval robust smoothing 



theory described in Ref. [28] for continuous-time uncertain 
systems satisfying a certain integral quadratic constraint. 

II. STANDARD QUANTUM LIMIT FOR 
ORNSTEIN-UHLENBECK PROCESS 

The standard quantum limit plays an important role as a 
benchmark for the quality of a measurement and is set by 
the minimum error in phase estimation that can be obtained 
using perfect heterodyne technique, or in other words, a 
non-adaptive filtering scheme. The case when the signal 
phase varies as a Weiner process has been considered in 
Ref. [11], where the minimum variance was obtained to be 
\JH I \\f2\ct\) . In this section, we shall consider the case when 
the signal phase varies as an OU process as in Ref. [23]. We 
deduce the minimum error covariance for the case of OU 
noise using the standard optimal filtering approach rather 
than the method used in Ref. [11]. 

The OU noise process under consideration is [23], [24]: 



4>{t) = -\<t>(t) + V^v(t), 



(i) 



where <fi(t) is the system phase to be estimated, A > is 
the mean reversion rate, k > is the inverse coherence time 
and v(t) is a zero-mean white Gaussian noise with unity 
amplitude. 

In our analysis here, we use the fact that the heterodyne 
scheme of measurement is, in principle, equivalent to, and 
incurs the same noise penalty as, dual-homodyne scheme 
[23], such as the schematic depicted in Fig.Q] We model the 
OU process as a signal at the input being phase-modulated 
using an electro-optic modulator (EOM) that is driven by an 
OU noise source. The modulated signal is then split using a 
50 — 50 beamsplitter into two arms each with a homodyne 
detector (HD1 and HD2, respectively, with the LO phase 
of HD1 7r/2 out of phase with that of HD2). The ratio of 
the output signals of the two arms goes to an arctan block, 
the output of which is fed into a low-pass filter (LPF). The 
filter for the case in Ref. [23] is sub-optimal, and we use 
here a Kalman filter, instead, to obtain the minimum error- 
covariance that determines the SQL for the OU process. 

The output signals of the two arms are: 



1 



Ji = -j= (2|a|sin< 



+ n 1 + n 2 ) , 



h = —/= (2|a| cos0 + ti 3 - n 4 ) , 
v2 

where ri\ and rt 3 are measurement noises of the two homo- 
dyne detectors, respectively, and n 2 and are the noises 
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Fig. 1. Block diagram of the dual-homodyne scheme for deducing the 
SQL for OU noise. 



arising from the vacuum entering the empty port of the input 
beamsplitter corresponding to the two arms, respectively. All 
these noises are assumed to be zero-mean white Gaussian. 
The output of the arctan block is: 

' 2\a\ sin tfi + n\ + n 2 



= arctan 



(2) 



2\a\ cos0 + ?t-3 — ri4 

A Taylor series expansion up to first-order terms of the 
right-hand side yields: 



2 a 2 a 



(3) 



The transfer function of the low-pass filter for the case of 
Ref. [23] is [24]: 



G(s) 



X 



(4) 



■& s + x 

We can determine the error covariance when that filter is 
used as follows. The system augmented with the filter may 
be represented by the state-space model: 

x = Ax Bw. 



(5) 



where 



and 
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From ([TJ, (f3]l and ©, we get: 



Process: 
Filter: 



<j> = -\<t> + 



' KV, 



2 a 



Thus, we have: 
A 



' -A 




and B = 


X ~X _ 
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2|a| 



n 2 . 



(6) 
(7) 








The steady-state state covariance matrix P5 is obtained 
by solving the Lyapunov equation: 

-t r 



AP S + P S A +BB 
where P5 is the symmetric matrix 

P s = £(xx T ) = 



0. 



(8) 



Pi P2 
P2 Ps 



Upon solving ®, we get 

P - K 
Pl - 2A' 

P2 = 



XK 



ft = 



2A(A + X )' 

X ( K 



1 



A(A + X ) 2|«|= 
The estimation error can be written as: 

e = <f> — 4> = [1 — 

which is mean zero since all of the quantities determining e 
are mean zero. 



The error covariance is then given as: 



a 2 = E(ee T ) = [1 - l]£(xx T ) 

1 



= [1-1 
Thus, we obtain: 

a 2 



Pl 
P 2 



Pi 
P 3 



-1 



1 

-1 
= Pi 

X 



2P 2 + P 3 



4 a 



(9) 



2(A + X ) 

Note that when A = 0, the above expression for the error 
reduces to that given by Eq. (3.8) from Ref. [11]. 

By contrast, the optimal filter is given by the Kalman filter, 
which may be determined in steady-state from the algebraic 
Riccati equation, which for the process given by (fl~|i and 
measurement given by (O is: 



-2AP -2\a\ 2 P 2 + k = 0. 



(10) 



The stabilising solution of the above equation for P is: 



P 



-A+ v /A 2 +2 K |q|- 
2\a\ 2 



(11) 



This, being the minimum error that can be obtained 
without feedback, determines the standard quantum limit 
for OU noise process. Note that when A — 0, we get 
P = (\/2\a\), as expected. 

The error-covariance given by (|9} is obtained when using 
the filter from Ref. [23], that is only optimal when assuming 
Wiener noise but not so for the more general OU noise, 
since it has only one variable \ that controls both the gain 
and the corner frequency of the filter. On the contrary, in 
deriving the error covariance given by (ITO of the Kalman 
filter, no such assumptions were made, so that the gain 
and the corner frequency of the filter were allowed to be 
two independent values, thereby yielding lower mean-square 
error in the estimate than in the former case. 

III. ADAPTIVE PHASE ESTIMATION USING 
SMOOTHER FROM REF. [23] 

In this section, we consider the (offline) estimator from 
Ref. [23], which is essentially a combination of two filters, 
one forward-time and another reverse-time. 

A. Forward Filter 

The forward filter has the following form [24]: 

Q-(t)= X f e(()e^- t U( = X (0(t)*e~x- t ), (12) 

J — oc 

where the value of \ is Xopt = 2|a|y / K, which is optimal 
in the limit A —> 0. Also, 9(t) comprise the measurements 
given by: 

m = m + -}-Mt), (B) 

2\a\ 

where w(t) is a zero-mean white Gaussian noise with unity 
amplitude. 



The transfer function of this filter is: 

e_( s ) x 



9(s) s + x 
Thus, the forward-time process and filter equations are: 



(14) 



Process: <p — ~M + V^u, 
Filter: 6>_ = -Xopt®- + Xo P *<5 



Xopt r 

2ld'' 



(15) 
(16) 



Let the steady-state state covariance matrix for the forward 
system be: 

S Mf 
M f Nf 



fs = 



(17) 



Note that 
Pfs =E 



Thus, we have 

Y. = E\<b\ M 



[ 4> e_ ] \=E 



06 _ 

e 2 _ 



// = £[0e_], N f = E[Q±}. 

Upon solving the Lyapunov equation of the form (O for 
the forward system, we get: 



2A' 



M f = 



XoptK 



2A(A + xopt) ' 

iXopt\a\ 2 K + XoptX 2 + xlpA 



8\a\ 2 X(\ + Xopt) 



Thus, we obtain: 



Xopt (A + 2xo P t) _ \/k(\ + 4|q|yg) 
8|a| 2 (A + Xo P t) ~ 4|a|(A + 2|a|v^)' 



(18) 
(19) 

(20) 



(21) 



One can verify that this expression for the error covariance 
agrees with a 2 _ of Eq. (10) from Ref. [23] for the optimal 
case of x- 

B. Backward Filter 

The backward filter has the following form: 

/>oo 

e+(t) = x / e(()e-x-«-Vdc 



(22) 



Let C = T - C and r = T - t. Then, d( = -d(. When 
C = t, ( = T — t = r and when C = oo, ( — T — oo = — oo. 
Thus, we get: 



e+(r) = -x 

6+(r)=x 



^)e-x.(r-C)^ = x ^ (r) 

where 6(() = 6{T - (). 

Thus, we get in the Laplace domain: 

e+(.s) _ x 
e(s) s + x' 

where s is the Laplace variable. 



(23) 



Thus, we obtain: 

6 + (r) 



where again the value of x is Xopt = 2|a|y / K, which is 
optimal in the limit A — > 0. 

When our model (fIll. dT3T l. which is driven by Gaussian 
white noise, has reached steady state, the output process will 
be a stationary Gaussian random process, which is described 
purely by its auto-correlation function. If we consider this 
output process in reverse time, this will also be a stationary 
random process with the same auto-correlation function. This 
follows from the definition of the auto-correlation function. 
Hence, the statistics of the reversed time output process are 
the same as the statistics of the forward time output process. 
Thus, the reversed time output process can be regarded as 
being generated by the same (and not time reversed) process 
that generated the forward time process, i.e. 

<Hr) = -A^(t) + Vkv(t), 

6(t) = <j>{r) + ^-Mt). 

2\a\ 

Now, we would use the Lyapunov method to deduce the 
state covariance matrix and error-covariance of the backward 
filter, the process and filter equations of which are: 



Process: </> = -X<f> + y/nv, 
Filter: 8+ = -Xo P t®+ + Xo P t<, 



Xopt 

2\a\' 



(24) 
(25) 



Let the steady-state state covariance matrix for the back- 
ward system be: 



bs = 



£ M b 
M b N b 



(26) 



where we have 

Z = E{<j) 2 }, M b = E[<j)Q + } 



m = Em 



Upon solving the Lyapunov equation of the form ^ for 
the information filter, we get: 

K 

2A' 

XoptK 



£ 

Mb = 



N h 



2A(A + xopt) ' 

4Xo P tH 2 K + Xopt A 2 + xlpA 



(27) 
(28) 



|a| 2 A(A + Xopt) 



Thus, we obtain: 



Xopt (A + 2xo P t) _ y^A + A\a\y/K) 
8|a| 2 (A + xopt) ~4H(A + 2|a|v^)' 



(29) 



(30) 



One can verify that this expression for the error covariance 
agrees with er + of Eq. (10) from Ref. [23] for the optimal 
case of x- 



C. Smoothed Error Covariance 

Suppose we have two unbiased estimates of some state x. 
We call these £\ and x 2 . We form a new estimate i as a 
linear combination of x\ and i 2 [29]: 

x = fcixi + k 2 X2, (31) 

where for the new estimate to be unbiased, 

fci + fc 2 = 1. 

The mean-square error for x is then 

E[{x-xf] =£{[x-fciii-(l-fci)i 2 ] 2 }, (32) 

or, 

E[e 2 ] =£{[ fcl ( ei - e2 ) + e2 ] 2 } 

= klE[el] + (1 - k-ifElel] + 2fci(l - ki)E[ ei e 2 }, 

where e, e\, and e 2 are the errors in x, X\, and £2, 
respectively. Eq. d33l may now be differentiated with respect 
to ki and set equal to to find the optimal k%. Thus, we 
obtain: 



E[el] 



E[eie 2 ] 



E[e\] + E{4] - 2E[ ei e 2 } ' 
Substituting for k\ in Q3t . we get: 



, = E[e\]E[el\ - (Eje^}) 2 
1 1 E[e 2 ]+E[e 2 2 ]-2E[e 1 e2\ 



(34) 



(35) 



This relation can be used to obtain the smoothed error 
covariance, given the forward and backward systems. In our 

case, e x = </> - 6_, e 2 = cj) - 6+, E[e\] = a}, E[e 2 2 ] = a 2 
and it remains to evaluate _E[eie 2 ] — E[(<p— 8_)(0 — 6+)]. 



E[ ei e 2 ] = E {[ 1 - 1 ] 
= [ 1 -1]E 
= \ 1 -1 



1 

-1 



06+ 





1 




1 — 1 



s M b 1 r 1 

Mf a^P \ [ -1 
= 11- Mf - M b + aS/8, 
where a = Af/S" 1 and /3 = E^Mb [16]. Thus, we get: 



( i 2 b = J B[(0-e_)(0-e + )] 



kA 



2(A + X opt) 2 



(36) 



Note that this agrees with Eq. (11) from Ref. [23] for the 
optimal case of x- 

Upon substituting appropriately in (l35l . we thus get the 
following as the smoothed error covariance: 



o 2 = 



Vk(A 2 +8|a|AV^ + 8|a| 2 K) 
8|a|(A + 2|a|^) 2 ' 



(37) 



One can verify that this expression for the error covariance 
agrees with a 2 of Eq. (12) from Ref. [23] for the optimal 
case of x- 



IV. ADAPTIVE PHASE ESTIMATION USING A 
RAUCH-TUNG-STRIEBEL SMOOTHER 

In this section we design the optimal RTS smoother for the 
adaptive system of Ref. [23] and analyse the error-covariance 
of the same to show that it is equivalent to an optimal two- 
filter smoother, before we compare it with the estimator used 
in Ref. [23] and the Kalman filter used in Ref. [24] in the 
next section. 

A. Forward Filter 

The forward filter is the same as the Kalman filter from 
Ref. [24]. The steady-state Riccati equation is: 



-2XP f - A\a\ 2 P] + k = 



(38) 



The stabilising solution of the above equation for Pf is: 

(39) 



^ = 4i^(-W4«M 2 + A*) 



The Kalman gain is: 



K f = -A + y/4n\a\ 2 + \ 2 . 
The forward filter equation is: 



= -(\ + K f )4> f +K f 0+-±-w. 

2 a 



(40) 



(41) 



B. Rauch-Tung-Striebel Smoother 

The steady-state Riccati equation for the RTS smoother is 
[22]: 



-2XP + 2kPJ 1 P - k = Q. 



(42) 



Upon substituting for Pf and solving for P, we get: 



2y/4 K \a\ 2 + X 2 



(43) 



Note that when A = 0, we get P — ^/k/A\q\, as desired 
(see Ref. [23]). 

The smoother gain is obtained to be [22]: 



F = 



4|a| 2 K 



-A+ ^An\a\ 2 + \ 2 
The equation for the smoothed estimate is [22]: 



(44) 



<f> = (-A + F)4> - F<f>f. 



(45) 



C. Backward Filter 

The RTS smoother, as obtained in the previous section, 
abstracts away the backward filter, the error-covariance and 
the filter equation of which are explicitly derived here for 
reference later. 

The steady-state Riccati equation for the backward filter 
is [22]: 



2XP b - A\a\ 2 P b 2 + k = 0. 



(46) 



p b = 7T12 f A + V^kH 2 + X' 



The stabilising solution of the above equation for Pi, is: 

(47) 



(48) 



Thus, the Kalman gain for the backward filter is: 



K b = A + ^An\a\ 2 + \ 2 . 
The backward filter equation is: 



<t>b = (A - K b )(j) b + K b (f> + -rKw. 

2 a 



(49) 



Remark.: One can verify that the error-covariances of 
the forward and the backward Kalman filters, obtained using 
the Lyapunov method used in section [TTT1 would be the same 
as in ( f39l > and d47b , respectively. In addition, referring to 
section ITlI-Cl one can show: 



E[(<j>-<l>f)(<j>-<k)]=0, 



(50) 



which implies that the forward and the backward estimates of 
the optimal smoother are independent. The error covariance 
of the RTS smoother thus obtained from ((35} would agree 
with ||43}. 

V. COMPARISON BETWEEN RTS SMOOTHER, 
KALMAN FILTER USED IN REF. [24], AND FILTER 
AND SMOOTHER USED IN REF. [23] 

Fig. |2] shows the plot of the mean-square error against the 
parameter A for the four cases, viz. RTS smoother, Kalman 
filter used in Ref. [24], filter used in Ref. [23] and smoother 
used in [23], as compared to the SQL for the noise process 
modulating the signal phase being estimated. The nominal 
experimental values used in the adaptive experiment in Ref. 
[23] were used for the other parameters in obtaining these 
graphs. It is clear that smoothing offers improvement over 
filtering alone in the accuracy of the estimate. For lower 
values of A, both the filter and the smoother from Ref. [23] 
approximate well the Kalman filter and RTS smoother, re- 
spectively. However, with increasing A, the optimal filter and 
smoother improve significantly and beat the SQL throughout, 
unlike the sub-optimal ones used in Ref. [23]. In summary, 
it is evident that the RTS smoother is the best scheme. The 
red vertical line denotes the value of A used in the adaptive 
experiment in Ref. [23]. 





— Filter from Ref. [23] 




— Kalman Filter from Ref. [24] 




— Smoother from Ref. [23] 

— RTS Smoother 
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Fig. 2. Comparison of the error covariance between the RTS smoother, 
Kalman filter used in Ref. [24], and filter and smoother used in Ref. [23]. 



VI. ADAPTIVE PHASE ESTIMATION USING A 
ROBUST FIXED-INTERVAL SMOOTHER 

We shall use the technique laid down in Ref. [28] to 
build a robust smoother for our continuous-time uncertain 
system, satisfying a certain Integral Quadratic Constraint as 
in Eq. (2.4) in Ref. [28]. The uncertainty is introduced in the 
parameter A as follows: 

A -> A - ^AA, 

where A is an uncertain parameter satisfying |A| < 1, and 
< fi < 1 determines the level of uncertainty in the model. 
Eq. (2.5) in Ref. [28], then, takes the form: 



Process: 
Measurement: 



= -\<f> + BxAK(t> + Btv, (51) 
1 



2ld 



(52) 



where B x = ^fn and K = [iXj ' %J~k. 

The uncertainty output of Eq. (2.1) in Ref. [28] for our 
system is: 

z = —=<p. 



(53) 



For the purpose of the Integral Quadratic Constraint (IQC) 
satisfied by our system, we will have Xq = 0, since no a- 
priori information exists about the initial condition of the 
state in our case. Also, Q = 1 for the uncertainty matrix 
A to satisfy the required bound. Also, d = 1, since the 
amplitudes of the white noise processes v and w have been 
assumed to be unity. The IQC of Eq. (2.4) in Ref. [28], thus, 
takes the following form in our case: 

/ (w 2 + -^ T ri 2 R)dt<l + [ \\z\\ 2 dt, (54) 
Jo 4|a| 2 J 

where w — AKcp + v and v — w are the uncertainty inputs. 

Thus, we would have R = 4|a| 2 in our case. 

The steady-state forward Riccati equation, as obtained 
from Eq. (5.1) in Ref. [28] for our case, is: 



2XX + kX 1 



4|a| 



0. 



(55) 



The stabilising solution of the above equation for X is: 

x = A+VA2-^ A 2 +4H 2^ 



The steady-state backward Riccati equation, as obtained 
from Eq. (5.2) in Ref. [28] for our case, is: 



- 2XY - kY z - 



4 \aV = 0. 



(57) 



The stabilising solution of the above equation for Y is: 

y = -X+^X 2 -^X 2 +4\a\ 2 K _ (5g) 



Next, Eq. (5.3) in Ref. [28] for our case yields: 

V = -(V x2 - M 2 ^ 2 + 4|a| 2 K)77 + 4|ct| 2 + 2\a\w. (59) 

Likewise, Eq. (5.4) in Ref. [28] for reverse-time yields: 

i = -W>? - M 2 A 2 + 4|a| 2 KK + 4|a| 2 + 2\a\w. (60) 

The forward filter is, then, simply the centre of the ellipse 
of Eq. (3.3) in Ref. [28]: 



Likewise, the backward filter is: 

The forward filter differential equation is, thus: 



= -L<f>f 



X + L ( 



2\a\n 

1 

X + L 



and the backward filter differential equation is: 



-L(j)b 



A\a\ 2 n 

-x + l' 



2\a\K 
-X + L 



(61) 



(62) 



(63) 



(64) 



where L = ^/X 2 — /i 2 A 2 + 4|a| 2 K. 

The robust smoother for the uncertain system would, then, 
be the centre of the ellipse of Eq. (5.5) in Ref. [28]: 



n - 



(65) 



X + Y 

One can verify that for [i = 0, d63l and d64b reduce to 
(PTlT i and (1491 . respectively. 

VII. COMPARISON BETWEEN ROBUST AND RTS 
SMOOTHERS FOR THE UNCERTAIN SYSTEM 

The error-covariances of the robust smoother and the 
RTS smoother for the uncertain system may be computed 
using the Lyapunov technique employed in section [in] as a 
function of A, for the nominal experimental values of all 
the parameters and a given value of /i. Fig. [3] shows the 
comparison of the two for the case of 50% uncertainty, below 
which the robust smoother is not quite significantly superior 
in performance as compared to the RTS smoother. 

Figs. @] and [5] show the comparison for /.i = 0.8 and p, = 
0.9, respectively. Clearly, the robust smoother performs much 
better than the RTS smoother as A approaches 1 for all levels 
of uncertainty in A. 



VIII. CONCLUSION 

This paper extends the optimal and robust filtering theory 
of Ref. [24], as applied to adaptive continuous homodyne 
phase estimation of a coherent state of light, to include 
optimal RTS and robust fixed-interval smoothing rather than 
filtering alone. In particular, it presents an insightful analysis 
of the relative performance of these various schemes with 
respect to the standard quantum limit. These theoretical 
results are to be demonstrated experimentally as part of 
further work. It would be interesting to extend these results 
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Fig. 3. Comparison of error covariance as a function of A for fi = 0.5. 
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Fig. 4. Comparison of error covariance as a function of A for fi = 0.8. 
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Fig. 5. Comparison of error covariance as a function of A for fj, = 0.9. 



for the case of squeezed states of light or other complex noise 
processes. Robustness to uncertainties in other parameters 
such as the photon flux or the noise power may also be 
explored. 
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